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ABSTRACT 


The  errors  generated  by  the  recovery  of  data  acquired  by  a  voltage - 
to-frequency  (V/F)  data  acquisition  system  are  theoretically  analyzed, 
using  a  digital  computer.  The  input  waveforms  considered  are  the 
sinusoid,  the  step  function,  and  the  ramp  function.  In  each  case,  the 
well  known  "sampling  theorem"  of  C.  E.  Shannon,  viz,  two  samples 
per  period  of  the  smallest  period  present  in  the  signal  are  sufficient 
for  recovery  of  data  with  perfect  fidelity,  will  be  shown  to  be  inapplicable 
as  a  criterion  for  accurate  and  reliable  recovery  of  data  acquired  by  a 
V/F  system.  The  errors  are  of  two  categories:  One  is  the  error 
from  linear  interpolation  recovery,  and  the  other  is  from  the  fact  that 
data  points  do  not  always  coincide  with  the  input.  The  sinusoid  analysis 
reveals  that  a  sampling  rate  of  11.  4  summing  intervals  per  period  is 
required  to  assure  5  percent  or  less  error.  The  step  function  analysis 
shows  that  the  errors  are  independent  of  sampling  rate;  however,  its 
response  time  is  approximately  equal  to  the  summing  interval.  The 
analysis  of  the  ramp  function  indicated  that  the  ratio  of  ramp  rise  time 
to  summing  interval  must  be  five  or  greater  to  assure  5  percent  or  less 
error.  These  criteria  are  unique  to  each  case,  and  in  the  case  of  the 
sinusoid,  are  many  times  larger  than  the  criterion  of  the  "sampling 
theorem.  " 


iii 


AEDC-TR-68-T43 


CONTENTS 

Page 

ABSTRACT .  iii 

NOMENCLATURE .  vi 

I.  INTRODUCTION .  1 

II.  VOLTAGE -TO-FREQUENCY  CONVERSION  METHOD  .  .  2 

HI.  SINUSOIDAL  INPUT 

3.  1  Linear  Interpolation  Error  Analysis  of  Sinusoidal 

Input .  4 

3.  2  Data  Point  Error  Analysis  of  Sinusoidal  Input  ....  7 

3.3  Synchronization  Effects  on  Errors  of 

Sinusoidal  Input .  8 

3.4  Experimental  Correlation  of  Theoretical  Errors  .  .  .  11 

IV.  STEP  INPUT 

4.  1  Linear  Interpolation  Error  Analysis  of  Step  Input  .  .  15 

4.  2  Data  Point  Error  Analysis  of  Step  Input .  16 

4.  3  Time  Delay  Analysis  of  Step  Input . .  .  17 

V.  RAMP  INPUT 

5.  1  Linear  Interpolation  Error  Analysis  of  Ramp  Input  .  .  21 

5.  2  Data  Point  Error  Analysis  of  Ramp  Input .  21 

VI.  SUMMARY  AND  CONCLUSIONS 

6.  1  Summary .  22 

6.  2  Conclusions .  23 

REFERENCES  . .  24 

APPENDIX:  Computer  Programs 

I-I.  Computer  Program  to  Find  Emax  or  emax  for  a 
Sinusoidal  Input  as  m  Varies  from  1.  0  to  25.  0 

in  Increments  of  0.  1 .  27 

I-II.  Computer  Program  to  Find  the  Largest, 

Smallest,  and  Delta  Emax  or  emax  as  xjj 
Varies  from  0  to  180  deg  in  Increments  of 
1  deg  and  m  Varies  from  1  to  25  in  Increments 

of  0.  5 .  29 

I -HI.  Computer  Program  to  Compare  Theoretical 

and  Actual  emax  or  Emax  for  Sinusoidal  Input  .  .  31 

I-IV.  Computer  Program  to  Find  emax  or  Emax  for 
a  Constant  K  as  /3  Varies  from  0.  1  to 
12.  0  in  Increments  of  0. 1 .  32 


v 


AEDC.TR.68. 143 


ILLUSTRATIONS 


Figure  Page 

1.  Block  Diagram  of  LRF  V/F  Acquisition  and 

Recovery  System .  2 

2.  Example  Illustrating  the  Summing  Interval  Method 

of  Data  Representation .  3 

3.  Plot  of  m  versus  emaX  f°r  Sinusoidal  Input .  5 

4.  Input- Output  Waveforms  for  Sample  Problem  ....  6 

5.  Plot  of  m  versus  Emax  for  Sinusoidal  Input .  7 

6.  Input-Output  Waveforms  for  Step  Input .  14 

7.  Plot  of  K  versus  emax  for  Step  Input .  16 

8.  Plot  of  K  versus  Emax  for  Step  Input .  17 

9.  Plot  of  K  versus  Tp  for  Step  Input .  18 

10.  Input-Output  Waveforms  for  Ramp  Input .  20 

11.  Plot  of  emax  versus  (3  for  Ramp  Input  with 

K  Constant .  21 

12.  Plot  of  versus  8  for  Ramp  Input  with 

KConstaS3. .  22 


TABLES 


I.  Emax  Solutions  of  Synchronization  Effects  on  Sinusoidal 
Inputs . 


II.  emax  Solutions  of  Synchronization  Effects  on  Sinusoidal 
Inputs . 


III. 


emax  an<*  Emax 
Sinusoidal  Input 


Theoretical  and  Actual  Solutions  for 


9. 

10 

13 


NOMENCLATURE 


E  Error  at  a  data  point 

e  Error  at  any  point  using  linear  interpolation 


AEDC-T  R. 68. 143 


f 

Frequency,  Hz 

Int  (  ) 

Integer  part  of 

K 

Synchronization  parameter  for  step  and  ramp  input 

Kv 

Gain  of  V/F  converter  (pulses/ sec -v) 

m 

Summing  intervals  per  period  of  sinusoid 

N  . 

Output  counts  of  V/F  converter  in  r  seconds 

n 

Number  of  arbitrary  summing  interval 

T 

Period  of  sinusoidal  input,  sec 

td 

Time  delay,  sec 

t 

Time 

V 

Voltage 

a 

Argument  of  Eqs.  (11)  and  (21) 

13 

Ratio  of  ramp  rise  time  to  summing  interval 

A<  ) 

Change  in 

Ui 

Step  function 

fj.  2 

Ramp  function 

T 

Length  of  summing  interval,  sec 

0 

Synchronization  parameter  for  sinusoidal  input 

SUBSCRIPTS 

i 

Input 

max 

Maximum  value 

min 

Minimum  value 

n 

Value  in  summing  interval  number  n 

o 

Output 

vii 


A  EDC-T  R-68-143 


SECTION  I 
INTRODUCTION 


Data  acquisition  relies  heavily  on  the  process  of  sampling  and  repre¬ 
senting  a  continuous  time  function  in  terms  of  discrete  time  samples. 
Numerous  methods  and  systems  to  accomplish  this  are  in  existence  and 
have  their  own  uniquenesses.  However,  all  of  these  methods  and  sys¬ 
tems  have  a  common  basis  that  justifies  the  representation  of  a  continu¬ 
ous  variable  by  discrete  samples,  and  that  is  the  sampling  theorem. 

The  sampling  theorem,  though  known  and  used  much  earlier,  is 
widely  attributed  to  C.  E.  Shannon  (Ref.  1).  It  may  be  stated  as  follows: 

"2/T  samples  per  second  suffice  to  represent  perfectly  and 
permit  perfect  recovery  of  a  time  function  provided  that  the 
time  function  contains  no  periodic  components  of  period  less 
than  T  seconds. " 

However,  the  data  recovery  from  an  integrating  system  such  as  the 
voltage -to-frequency  data  acquisition  system  when  used  as  a  sampling 
system,  i.  e. ,  used  to  recover  the  actual  waveshape,  does  not  have  a 
direct  application  through  the  sampling  theorem. 

This  investigation  is  a  derivation  of  the  number  of  integrating 
intervals  required  to  accurately  reconstruct  an  input  using  linear  inter¬ 
polation  of  the  data  points  obtained  from  an  integrating  data  acquisition 
system.  Assuming  that  an  error  also  exists  at  the  data  points,  the 
analysis  will  include  both  a  determination  of  data  point  error  and  linear 
interpolation  error.  Error  will  be  the  absolute  value  of  the  difference 
of  the  input  and  the  representative  output  at  a  given  time  as  a  percentage 
of  the  maximum  magnitude  of  the  input.  This  is  a  modification  of  the 
definition  of  deviation  as  defined  in  Ref.  2,  page  209. 

Three  ideal  signals  will  be  treated  theoretically:  sinusoid,  step, 
and  ramp.  The  digital  computer  will  be  used  to  simulate  this  type  of 
data  recovery,  and  the  maximum  error  for  various  parameter  configura¬ 
tions  will  be  determined.  From  this,  the  limits  of  the  parameters  for 
accurate  (5  percent)  recovery  can  be  determined. 

Experimental  verification  of  the  error  analysis  for  the  sinusoid  was 
performed  using  the  VIDAR®  voltage -to-frequency  data  acquisition  sys¬ 
tem  in  the  Large  Rocket  Facility  at  AEDC. 
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SECTION  II 

VOLTAGE-TO-FREQUENCY  CONVERSION  METHOD 


The  method  of  data  acquisition  that  is  to  be  treated  in  this  study  is 
a  voltage -to -frequency  (V/F)  conversion.  The  block  diagram  in  Fig.  1 
depicts  the  setup  of  such  a  system. 


TAPE 

RECORDER 


Fig.  1  Block  Diagram  of  LRF  V/F  Acquisition  and  Recovery  System 


The  V/F  converter  develops  output  pulses  at  a  rate  precisely 
proportional  to  the  input  voltage.  The  output  pulses  are  recorded  on 
magnetic  tape  for  later  data  reduction.  The  recorded  signals  are 
recovered  by  means  of  a  translator,  which  works  like  a  gated  frequency 
counter  that  counts  the  pulses  during  a  time  interval,  r.  The  counts 
contained  in  a  summing  interval  are  converted  to  a  representative  volt¬ 
age,  VTn,  by  a  digital  computer  using  the  calibration  of  the  V/F  con¬ 
verter  and  summing  interval  of  the  translator. 

To  illustrate:  Let  the  V/F  converter  be  adjusted  so  that  the  output 
is  N0  counts  in  t  seconds  with  the  input  shorted.  Then  apply  a  calibra¬ 
tion  voltage  of  Vc  volts  to  the  input  of  the  V/F  converter;  let  Nc  be  the 
counts  in  r  seconds  of  the  output.  For  a  number  of  counts  in  the  nth 
interval  of  t  for  an  arbitrary  input  of  V^(t)  defined  as  NTn,  the  following 
expression  defines  VTn: 


'V-  = 


Nrn-N. 


V, 


(1) 


But,  because  the  V/F  data  recovery  technique  is  an  integrating 
method, 

/(n  +  Dr 

Vi<t)  dt  (2a) 


and 


I 


Nc.  -  N0  s  I  Kv  Vcdt  =  Kv  Vcr 


(2b) 
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where  Kv  is  the  gain  of  V/F  converter  with  units  of  pulses  per 
second-volt. 


Substituting  the  definitions  of  Eqs.  (2a)  and  (2b)  into  Eq.  (1),  it 
becomes: 


Vr„  =  t 


*1 


(n  +  1 )T 

”  Vj(t)  dt 


(3) 


Equation  (3)  will  be  used  to  theoretically  simulate  the  V/F  conversion 
method  of  data  acquisition  and  analyze  the  errors  resulting  from  taking 
the  value  of  vm  as  a  data  point  that  occurs  at  the  midpoint  of  the  slim¬ 
ming  interval  from  nr  to  (n+l)r  and  the  points  linearly  connected  to 
represent  the  input  waveform.  Figure  2  illustrates  this  type  of  data 
acquisition  and  waveform  reconstruction. 


Fig.  2  Example  Illustrating  the  Summing  Interval  Method  of  Data! Representation 
Since  the  value  Vrn  is  an  integral,,  its  value  and  the  value  of  Vj(t) 

2n+i 

at  t  =  Tn  =  — g — t,  the  time  at  the  midpoint  of  the  nth  interval,  are  not 

always  equal;  therefore,  the  error  analysis  also  investigates  the  data 
point  error.  To  distinguish  the  two,  the  lower  case  letter,  e,  desig¬ 
nates  the  error  from  linear  interpolation  and  the  upper  case  letter,  E, 
is  used  to  designate  the  error  at  data  points. 
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SECTION  III 
SINUSOIDAL  INPUT 


2  7Tt 

For  the  sinusoidal  analysis,  let  V^(t)  =  sin  which  is  normalized 
for  ease  of  application  and  does  not  affect  the  actual  error  analysis. 

Substituting  this  function  in  Eq.  (3),  the  expression  for  the  repre¬ 
sentative  output  becomes: 

/(n  +  l)r 

sin  J, 

'I 

E_  (n  +  Dr 

os^-l 

TJ„ 


--k-E 


2  77(11  +  l)r  2rrn  r 

cos  - - -  -  cos  — — 


(4) 


Using  the  identity  from  Ref.  3,  page  18  for  the  difference  of  two 
cosine  functions,  Eq.  (4)  becomes: 


\j  T  P  .  (2n  +■  1  )rrr  . 

Vr  »  —  sin -  •  sin  — 

"  trr  L  T  T,J 


(5) 


Letting  the  ratio  T/t,  which  is  the  number  of  summing  intervals  per 
period,  be  designated  with  the  letter  m,  the  expression  in  Eq.  (5)  becomes: 


„  m|  r  .  (2n  -  l)ff 

V-  =  —  sin  ■ — -  •  sin  - 

a  IT  L  m  ®J 


(6) 


3.1  LINEAR  INTERPOLATION  ERROR  ANALYSIS  OF  SINUSOIDAL  INPUT 

An  expression  for  e  consists  of  n  equations;  therefore,  the  deter¬ 
mination  of  a  maximum  value  is  a  tremendous  task  to  perform  manually. 
Consequently,  a  digital  computer  program  was  written  to  determine  the 
errors  for  a  sinusoidal  input  as  m  is  varied  from  1  to  25. 0  in  increments 
of  0. 1.  The  program  was  designed  to  search  the  number  of  cycles  that 
include  a  whole  number  of  summing  intervals.  The  reason  for  this 
approach  to  the  emaX  determination  is  that,  after  this  number  of  cycles, 
the  errors  become  repetitive.  By  incorporating  this  test,  the  emax 
determination  is  expedited.  The  program  is  listed  in  Table  I-I 
(Appendix). 
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The  solutions  for  emax  versus  m  are  plotted  in  Fig.  3.  The  plot 
depicts  an  irregular  relationship,  and  the  most  significant  fluctua¬ 
tions  are  in  the  range  of  m  =  1.0  to  3.0.  Therefore,  little  faith  can 
be  placed  in  the  application  of  the  basic  sampling  theorem  of  two  sum¬ 
ming  intervals  per  period  to  the  V/F  system.  Another  significant 
factor  is  that,  when  the  sampling  rate  is  close  to  frequency  components 
of  the  data  that  have  prominent  magnitude,  the  errors  greatly  reduce 
the  accuracy  of  the  data  recovery. 


As  an  example  of  the  use  of  Fig.  3,  consider  an  input  voltage  signal 
that  contains  a  50 -Hz  component.  What  error  would  be  introduced  in 
the  data  reduction  of  this  component  if  r  were  10  msec?  What  is  the  error 
for  r  =  2  msec?  For  r  =  10  msec. 


io3 

m  =  - 

.(50)  (10) 

From  Fig.  3,  emax  =  44  percent. 


2 


For  r  =  2  msec. 


m 


103 
(50)  (2) 


From  Fig.  3,  emax  =  6  percent. 
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Fig.  3  Plot  of  m  versus  emsx  for  Sinusoidal  Input 
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Figure  4  depicts  the  input -output  relationship  for  these  two  summing 
intervals  for  the  50 -Hz  component.  . 

In  order  that  emax  be  5  percent  or  less  for  a  sinusoidal  input,  m 
must  be  greater  than  11.4  and  at  m  =  2,  emax  =  46  percent.  This  indi¬ 
cates  a  large  error  at  the  rate  given  by  the  sampling  theorem. 


Fig.  4  Input-Output  Waveforms  for  Sample  Problem 


6 


AEDC-TR- 68-143 


3.2  DATA  POINT  ERROR  ANALYSIS  OF  SINUSOIDAL  INPUT 

The  error  analysis  for  Emax  is  much  simpler  than  for  emax.  The 
data  point  error  is  defined  by  this  relationship: 

E  =  100  | Vi  (rn)  -  VrJ  (7) 

since  |Vj(t)|max  is  equal  to  one. 

Combining  Eqs.  (6)  and  (7)  and  substituting  the  normalized  sinusoid 
with  the  definition  of  m  included,  the  expression  for  E  becomes: 

E  =  100  |sin  i2n-2_U"(l  _  £sin£-)|  (8) 

Even  though  the  error  is  expressible  in  one  equation,  many  calcu¬ 
lations  are  required  to  find  Emax>  Therefore,  the  program  used  to  find 
emax  a^so  ^as  9X1  ^max  option.  These  solutions  are  depicted  in  Fig.  5. 


Fig.  5  Plot  of  m  versus  Emox  for  Sinusoidal  Input 

This  plot  of  EmaX  versus  m  is  a  smoother  curve  than  the  one  for 
emax  versus  m.  The  minimum  value  of  m  for  all  data  points  to  be 
within  5  percent  is  5.  7,  and  at  the  basic  sampling  rate  given  by  the 
sampling  theorem  (two  summing  intervals /period),  the  data  point  is  in 
error  at  a  maximum  by  36+  percent. 
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3.3  SYNCHRONIZATION  EFFECTS  ON  ERRORS  OF SINSUOIDAL  INPUT 

The  analyses  of  the  sinusoid  for  emax  and  Emax  have  made  the 
assumption  that  the  sinusoid  and  the  n  =  0  interval' both  started  at  a 
time  when  the  value  of  the  sinusoid  was  zero  and  had  a  positive  slope. 
This  section  deals  with  the  effects  of  synchronization  upon  emax  and 

®max‘ 

The  input  signal  must  be  expressed  with  a  synchronization  param¬ 
eter,  ip.  The  input  then  is: 


V,(t)  =  sin^p-  + 

Substituting  this  expression  for  V^(t)  into  Eq.  (3): 


Using  the  definition  of  m,  Eq.  (9)  becomes: 

Vr  =  “  sin  pin. +  D ir  +  0-|  .  sb  51  (10) 

Since  many  solutions  of  Eq.  (10)  are  required  to  determine  the 
synchronization  effects,  a  program  was  written  to  vary  m  from  1  to 
25.  0  in  increments  of  0.  5  and  to  vary  ip  from  0  to  180  deg  in  increments 
of  1  deg.  The  solutions  for  the  Emax  option  are  in  Table  I  and  the  solu¬ 
tions  for  the  emax  option  are  in  Table  II.  The  program  is  listed  in 
Table  I-II  (Appendix). 
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TABLE  I 

Emox  SOLUTIONS  OF  SYNCHRONIZATION  EFFECTS  ON  SINUSOIDAL  INPUTS 


DATA  POINT  ERROR  ANALYSIS' 


LARGEST 

SMALLEST 

DELTA 

M 

PS  I 

‘  E 1  MAX ) 

PSI 

El  MAX) 

E(MAX) 

I.  00 

90.00 

100.00 

0.0 

0.00 

100.00 

1.50 

3C.00 

58.65 

0.0 

50.79. 

7.86 

?. .  OC 

0.0 

36.34 

90.00 

O.OO 

36.34 

2.5C 

18.00 

24.32 

72.00 

23.13 

1.19 

3.00 

30.00 

•  -  17.30 

.  Q.O  . 

14.98 

2.32 

3.5C 

90.00 

12.90 

0.0 

12.57 

0.32 

4.0C 

49.00 

9.97 

0.0 

7.05 

2.92 

4.50 

10.00 

7.93 

0.0 

7.81 

0.12 

5.00 

18.00 

6.45 

0.0 

6.  14 

0.32 

5.50 

41. CC 

5.35 

0.0 

5.30 

0.05 

6.00 

C.P 

4.5  1 

30.00 

3.90 

0.60 

6.50 

90.00 

3.85 

0.0 

3.82 

0.03 

7.00 

90.00 

3.32 

0.0 

3.24 

0.08 

7.50 

6  •  0  0 

2 .90 

60.00 

2.89 

0.02 

9.00 

22.00 

2.55 

90.00 

2.36 

0.19 

8.50 

16.00 

2.26 

180.00 

2.25 

0.01 

9.00 

10.00 

_2 .02 

0.0 

1.99 

0. 03 

9.5C 

52. CO 

1.81 

0.0 

1.81 

0.01 

10.00 

0.0 

1.64 

18.00 

1.56 

C.  08 

10.50 

73.00 

1.49 

120.00 

1.43 

0.00 

11.00 

8.00 

1.35 

180.00 

1.34 

0.01 

1  L.5C 

43.00 

1.24 

0.0 

1.24 

0.00 

12. CO 

15.00 

.1.14 

0.0 

1.10 

0.04 

12. 5C 

11.00 

1.05 

0.0 

i.05 

0.00 

13.00 

76.00 

0.97 

0.0 

0.96 

0.01 

13.50 

10.00 

0.90 

0.0 

0.90 

0.00 

14. OC 

1C3.00 

0.84 

90.00 

0.82 

0.02 

14.  50 

34.00 

6.73 

0.0 

0.78 

0.00 

..  .  15.00. 

6. CO 

0.73 

0.0 

0.73 

0.00 

15.50 

32.00 

0.68 

l’so.oo 

0.68 

0.00 

16.00 

11.00 

0.64 

45.00 

0.63 

0.01 

16.50 

101.00 

0.60 

C.O 

0.60 

0.00 

17.00 

16.00 

0.5  7 

0.0 

0.57 

0.00 

17.50 

85.00 

0.54 

108.00 

0.54 

0.00 

_._1.8,.QC.. 

0.51 

10.00 

0.50 

0.01 

18.50 

17.00 

6.48 

107.00 

6.48 

0.00 

19.00 

71.00 

0.46 

180.00 

0.45 

0.00 

19.50 

67.00 

0.43 

60.00 

0.43 

o.oo 

_ 2.0.00  . 

9.00  - 

0.41 

0..  0 

0.41 

0.01 

20.50 

11.  GO 

.  0.39 

180.00 

6.39 

0.00 

.21.00 

.  _  ,4. 00 

0.37 

180.00 

0.37 

0.00 

21.50 

23. CO 

0.36 

21.00 

0.36 

0.00 

22.00 

33.00 

0.34 

9C.00 

0.34 

C.00 

22.50 

2. CO 

0.32 

12.00 

0.32 

0.00 

23.00 

59.00 

0.31 

0.0 

0.31 

0.00 

23.50 

6.00 

0.30 

46.00 

0.30 

0.00 

_ -2Ju.no 

8,00 

0.29 

180.00 

0.28 

0.00 

24.50 

46.00 

6.27 

136.00 

6.27 

o.oo 

25.00 

4.00 

0.26 

108.00 

0.26 

0.00 
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TABLE  II 

omax  SOLUTIONS  OF  SYNCHRONIZATION  EFFECTS  ON  SINUSOIDAL  INPUTS 


M 

PM 

1.00 

1.50 

0  0.00 

2.00 

QO.ee 

2.50 

90.  CO 

3.00 

QC.OO 

3. 5C 

9C.C0 

4.CC 

90. CO 

4.50 

9  o.OO 

5.00 

162. CO 

5. 50 

9C.OO 

6.0C 

90.  OC 

6. 50 

90.00 

7.00 

90.00 

7.5C 

9C.CC 

8.00 

9C.CC 

8.50 

90.00 

0.00 

90.00 

0.50 

90.00 

10.00 

126.00 

10.50 

90.00 

11. GO 

90.00 

1  1.50 

9C.C0 

12. GO 

0.0 

12.50 

54.  GO 

1  3.00 

9C.00 

M.5C 

9C.00 

14.30 

90.  GO 

14.50 

90.00 

15.00 

18. CO 

15.50 

9G.00 

14.  OC 

0.0 

16.50 

90.00 

17.00 

99.00 

17.50 

90.00 

18.00 

9C.00 

18.50 

90.90 

IQ.  or 

90.00 

10.50 

90.00 

20.  CO 

0.0 

20.50 

9C.C0 

?  1.00 

90.  CO 

21.50 

9C.C0 

22.00 

90.00 

22.50 

90.00 

23.00 

90.00 

2  3.50 

90.00 

24.00 

0.0 

2  4.50 

90.00 

2  5.0C 

11.00 

L  ~,CST 


F( MAX) 

PSI 

ICC.  r.  c 

9.00 

120.6  V 

0.0 

ICO.  CO 

c.o 

76.61 

C .  9 

58.65 

C.O 

4  5.69 

C.O 

36.34 

4  5.00 

29.47 

C.O 

24.32 

3  6.00 

20.38 

164.00 

1  7.30 

O.C 

14.  86 

55. CO 

12.90 

0.0 

11.29 

C.C 

9.97 

68.00 

8.  86 

C.O 

7.93 

0.0 

7.13 

126.00 

6.45 

0.0 

5.66 

123.00 

5.35 

100.00 

4.90 

69.00 

4.51 

135.00 

«♦.  16 

108.00 

3.85 

127. CO 

3.57 

23.00 

3.32 

0.0 

3.10 

57.00 

2.90 

36.00 

2.72 

54.09 

2.55 

55.00 

2.49 

19.  00 

2.26 

0.0 

2.13 

108.00 

2.02 

104. or 

1.91 

0.0 

1.81 

79.  00 

1.72 

50. CO 

l  •  64 

117.00 

1.56 

20.00 

1.49 

33.  00 

1.42 

26.00 

1.35 

105.00 

1.29 

18C.00 

1.24 

23.00 

1.19 

0.0 

1.14 

43.00 

1.09 

0.0 

1.C5 

106.00 

SMALLEST 

DELTA 

E l MAX ) 

E  C  MAX ) 

99.95 

0.05 

113.77 

6.91 

44.65 

55.35 

73.76 

2.85 

52.15 

6.  50 

44.68 

1.02 

26.89 

9.45 

29. C3 

0.44 

23.18 

1.14 

20.15 

0.23 

15.06 

2.24 

14.74 

C.  12 

12.56 

0.33 

11.29 

C.  10 

9.24 

0.73 

8.79 

0.07 

7.79 

0.14 

7  .  C  3 

0.05 

5.  14 

0.  32 

5.81 

9.05 

5.27 

C .  08 

4.85 

0.05 

4.32 

0.19 

4.15 

0.01 

3.1) 

0.06 

3.53 

0.04 

3.19 

0.13 

3.C7 

9.0  3 

2.86 

C  •  94 

2.70 

C .  92 

2.47 

0.08 

2.33 

0.02 

2.21 

0.0  5 

2.12 

0.02 

1.97 

0.05 

1.90 

0.02 

1.79 

0.93 

1.71 

J .  0  2 

1.62 

0.92 

1.54 

0.02 

1.46 

0.02 

1.39 

C.  02 

1.31 

0.04 

1.26 

0.0  3 

1.20 

0.93 

1.15 

0.04 

1.09 

0.05 

1 .04 

C.  05 

1.C5 

C.00 
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These  solutions  depict  the  cases  in  which  the  variations  of  Emax 
or  emax  are  the  largest.  The  cases  are  the  integer  values  of  m  and 
m  with  a  decimal  part  of  0.  5.  All  other  values  of  m  exhibit  less  varia¬ 
tion  in  Emax  or  emax  as  ^  varies. 

From  Table  I, 

for  m  >  6,  AEmax  <  0.  19  percent 
for  m  >  11. 4,  AEmax  <  0.  04  percent 
for  m  >  14.0,  AEmax  <0.01  percent 

The  m  >  11.  4  values  yield  accurate  recovery. 

From  Table  II, 

for  m  >  8,  Aemax  <0.32  percent 
for  m  >  11.4,  Aemax<  0. 19  percent 
for  m  >  16,  Aemax  <  0.  05  percent 

This  clearly  indicates  that  synchronization  has  negligible  effect  on  the 
error  at  a  summing  interval-to-period  ratio  that  is  required  to  accu¬ 
rately  recover  the  data. 

3.4  EXPERIMENTAL  CORRELATION  OF  THEORETICAL  ERRORS 


A  verification,  by  actual  errors  encountered  in  the  data  recovery, 
of  the  theoretical  treatment  was  desired;  therefore,  an  experiment  to 
determine  actual  errors  was  formulated.  For  the  input  waveform,  a 
sinusoid  was  chosen.  There  were,  basically,  two  reasons  for  this 
choice:  (1)  a  sinusoid  is  easily  obtainable  from  a  sinusoidal  generator, 
and  (2)  the  frequency  or  period  of  cycle  is  a  very  important  factor  in  the 
error  of  a  signal. 

The  VIDAR  Data  Acquisition  System  was  set  up  with  a  sinusoidal 
generator  connected  to  the  input  of  a  channel  of  the  V/F  converter.  The 
frequency  output  was  recorded,  using  direct  electronics,  on  a  magnetic 
tape  recorder.  The  recorded  tape  was  translated  using  a  summing 
interval  of  2  msec,  and  the  raw  counts  were  printed  out. 

In  order  to  gee  the  negative  values  of  the  sinusoid,  the  VIDAR  was 
adjusted  so  that  for  zero  input,  the  output  counts  were  at  midscale. 


11 


AEDC-TR-68-143 


Then  a  known  voltage  of  1.  5  v  was  applied,  and  the  VIDAR.  spanned  to 
nearly  full  scale.  The  count  levels  of  these  adjustments  were: 

Ov  Input  =  25,  003  counts/sec 
1.  5v  Input  =  47,  500  counts/ sec 

The  signal  from  the  sinusoid  generator  was  set  at  20  Hz  and  2  v 
peak-to-peak  which  is  easily  accommodated  in  magnitude  by  the  chosen 
calibration  voltage. 

Since  the  theoretical  error  analysis  used  m  as  the  variable,  m  had 
to  be  determined  for  this  case  of  the  translation  summing  interval;  m 
was  defined  as  T/r  or  1/fr;  giving  m  =  1/20  (0.  002)  =  25. 

Synchronization  of  output  and  input  for  actual  error  analysis  proved 
to  be  a  difficult  problem;  therefore,  since  for  an  m  of  25,  emax  =  1.05 
percent  and  Emax  =  0.  26  percent,  this  linear  interpolated  curve  was 
used  as  the  reference  or  standard.  Small  m's  were  obtained  by  taking 
two  or  more  adjacent  summing  intervals  at  a  time.  A  computer  program, 
depicted  in  Table  I-III  (Appendix),  was  written  to  do  the  task  of  summing 
up  adjacent  summing  intervals.  The  program  limits  the  number  of 
adjacent  summing  intervals  to  13  since  this  gives  an  m  of  1.92,  which 
is  below  the  level  set  forth  by  the  sampling  theorem.  Table  III  is  the 
solution  for  this  program  which  has  the  option  of  determining  either 

emax  or  ^max- 

As  shown  in  Table  III,  the  maximum  errors  determined  by  the  theo¬ 
retical  approach  were  confirmed  to  be  reasonably  accurate;  therefore, 
the  general  rule  of  thumb  of  11.4  summing  intervals  per  period  for 
emax  —  5  percent  is  a  valid  one.  This  also  indicates  that  the  theoretical 
simulation  is  a  reasonably  accurate  error  determination  technique. 


SECTION  IV 
STEP  INPUT 


In  the  theoretical  analysis  of  errors  of  the  V/F  converter  data 
acquisition  system  in  following  a  step  input,  the  time  reference  of  t  =  0 
is  set  at  the  start  of  the  summing  interval  in  which  the  step  function 
occurs. 
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TABLE  III 

®mox  and  Emax  THEORETICAL  AND  ACTUAL  SOLUTIONS  FOR  SINUSOIDAL  INPUT 


ACTUAL 

THEORY 

ACTUAL 

THEORY 

11 

EC MAX) 

EC MAX) 

* 

M 

EC MAX) 

EC MAX) 

1.92 

103.31 

103.03 

2.03 

34.32 

95.63 

2.27 

35.  40 

30.51 

2.50 

7  5.63 

73.70 

2.73 

64.5  0 

05.31 

3.13 

54,  01 

54.39 

3.57., 

'  43.32 

44.03 

**  .  * 

4.17 

33.70 

33.75 

5.00 

23.33 

23.13 

6.25 

16.34 

15.93 

3.  33 

10.47 

9.19 

12.50 

5.35 

4.15 

DATA  POI NT  ERROR 

ANALYSES'  -V 

.  ACTUAL- 

theory-  ,  <  •  £<=; 

*  c 

ACTUAL  i. 

—THEORY  - 

M 

ECMAX) 

ECMAX) '  •' 

•  .-.r-  M'  ' 

ECMAX)  ' 

-ECMAX) 

1.92 

33.67 

33.33 

2. OS 

32.33 

33.75 

2.  27 

23.57 

23.33 

2.  50 

22.35 

23. 13 

2.73  . 

20.31 

19.96 . 

•  '  :3 . 1 3  < 

.15 . 35  * 

15. .93  . 

3.57 

13.47 

12.33-  '  : 

.  4.17  > 

---'8.57 

,  ■  9.1-9;' 

5. 00,! 

-  3  'i  0  O'  - 

:s  -  C .  14  ***••* 

■",J&  -uv25<-  1 

K'v-4.  29- 

.  -  -4  ,- 15-': 

3.33 

4.  76 

2.35 

12.50 

0.01 

1.05 

The  input  waveform,  a  normalized  step  function,  is  denoted  by  the 
symbol : 


A,  (t  -  a)  =  1  <  a 

)1  t  >  a 

(ID 

For  this  analysis,  the  input  is: 

Vj{t)  =  ^  (t  -  Kr) 

(12) 

where  0  <  K  <  1  is  the  shift  of  the  step  with  respect  to  the  summing 
interval  reference. 


Using  Eqs.  (3)  and  (12)  to  find  the  output,  the  particular  interval 
for  n  =  0  is  the  only  one  of  interest  since  in  all  intervals  for  n  <  0, 

VTn  =  °>  and  in  all  intervals  for  n  >  0,  VTn  =  1. 

The  interval  for  n  =  0  has  an  output  value  of: 

VrD  =  if  ix,  (t  -  Kr)  dt 

■  if*  -  U’ -K) 

=  1  -  K  (13) 
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The  input -output  relationship  is  depicted  in  Fig.  6. 


OiKil/2 


1/2  <  K  i  I 

Fig.  6  Input-Output  Waveforms  for  Step  Input 
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4.1  LINEAR  INTERPOLATION  ERROR  ANALYSIS  OF  STEP  INPUT 


As  shown  in  Fig.  6,  the  maximum  error  occurs  at  t  =  Kt.  Point  a 
is  used  to  designate  this  intersection  of  the  input  and  output  waveforms 
at  t  =  Kt.  The  output  line  from  t  =  -  ^  to  —  is  designated  VQ.  (t),  and 

u  u  1 


T  3  X 

the  output  line  from  t  =  ^-to  ^  is  designated  V0„  (t).  For  0  <  K  <  1/2, 


point  a  occurs  on  the  output  line  VQ1  (t).  The  equation  for  this  line, 

1  X 

which  has  a  slope  of  —  (1-K)  and  passes  through  the  point  (-— ,  0),  is: 

X  b 


Since  (Vj  (t)|  max 


V0  (t)  =  -  (1  -  K)  (t  +  9 

IT  2 


emax  100  ^  o(  (Kt) 

=  100  (L  -  K)  (K  +  lA) 

-  SO  (l  -  K  -  2K2)  (14) 

for  0  <  K  <  1/2. 

For  1/2  <  K  <  1,  point  a  occurs  on  the  output  line  V02  (t).  The 
equation  for  this  line,  which  has  a  slope  of  K/r  and  passes  through  the 

point  (|^,  1),  is: 


Likewise, 


for  1/2  <  K  <  1. 


Vo,  (t)  -  (t  -  2p  +  1 

em  ax  =  100  []  -  V0j  (Kr)] 
=  -100K  (K  -  3/2) 

=  150K  -  100K2 
=  50  (3K  -  2K2) 


(15) 


Figure  7  is  a  plot  of  emax  versus  K  using  Eqs.  (14)  and  (15).  As 
can  be  seen,  emax  varies  from  the  lowest  value  of  50  percent  at  K  =  0, 
1/2,  and  1  to  a  highest  value  of  56.  25  percent  at  K  =  1/4  and  3/4. 
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Fig.  7  Plot  of  K  versus  emax  for  Step  Input 


4.2  DATA  POINT  ERROR  ANALYSIS  OF  STEP  INPUT 


For  this  input  signal,  all  data  point  errors  are  0  percent  with 
the  exception  of  the  one  for  n  =  0.  This  data  point  of  interest  occurs 
at  time  t  =  t0  or  0.  5. 


Using  Eq.  (13)  and  the  fact  that  jvi(t)  |  max  =  1,  the  Emax  for  step 
input  is: 


Emax  =  100  |  V j  (r0)  -  (]  -  K)| 


(16) 


However,  Vi(r0)  is  dependent  on  the  value  of  K  in  the  following 


manner: 


Combining  Eqs.  (16)  and  (17), 

(]  00K 

Emax  =  boo  (i 


0  <  K  <  0.5 

0.5  <  K  <  1.0  (17) 

the  maximum  error  becomes: 

o  <  K  <  0.5 
.  K)  0.5  <  K  <  1.0 
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A  plot  of  Emax  versus  K  is  shown  in  Fig.  8.  As  was  the  case  in  the 
total  error,  the  maximum  data  point  error  is  independent  of  the  summing 
interval  length  and  is  a  function  of  K  alone. 


Fig.  8  Plot  of  K  versus  Emox  for  Step  Input 


4.3  TIME  DELAY  ANALYSIS  OF  STEP  INPUT 

The  determination  of  time  delay,  i.  e. ,  time  for  the  response  of 
the  linear  interpolated  output  to  be  within  5  percent  error,  was  desired. 
The  notation  Tj)  is  used  for  this  parameter. 

This  determination  must,  like  the  error  analysis,  be  divided  into 
two  parts.  Likewise,  the  first  part  uses  the  output,  V01  (t),  in  its  deter¬ 
mination.  The  ranges  of  K,  however,  are  not  the  same.  In  using  VG1  (t), 
the  following  inequality  must  be  satisfied: 

1  -  V'o/f)  <  °.°5  <18> 

Since  V0^  (^-)  =  1-K,  Eq.  (18)  defines  the  limits  of  K  to  be: 

0  <  K  <  0.05 
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For  a  time  delay  of  Tjj,  the  actual  time  t  =  Tjj  +  Kt,  and  from  the 
equation  for  VQ^  (t), 

v0l  (To  t  Kr)  =  0.95  =  j  (1  -  K)  (Td  +  Kr  +  y) 


Solving  for  Tp, 


for  0  <K  <0.05. 


To 


0.45  -  0.5K  +  K2 
1  -  K 


(19) 


The  second  part,  0.05  <K  <  l,  uses  VQ2  (t)  in  the  determination  of 
Td,  with  t  =  Td  +Kt. 

Vq2  (Td  +  Kr)  =  0.95  =  (Tq  +  Kr  -  1.5r)  +  1 


Solving  for  Td, 

Tn  =  ^°-05  +  1.5K  -  K2 
Ka 

for  0.05  <K<  1. 


(20) 


Using  Eqs.  (19)  and  (20),  a  plot  of  Td  versus  K  is  depicted  in 
Fig.  9.  As  shown  in  this  plot,  Td  varies  from  a  minimum  of  0. 449r 
to  a  maximum  of  1.  053t. 


Fig.  9  Plot  of  K  versus  Tq  for  Step  Input 
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The  step  response  will  always  appear  to  be  a  ramp.  For  K  =  1/2, 
the  rise  time  of  the  ramp  is  2t  sec,  and  for  K  =  0  or  1,  the  rise  time 
will  be  r  sec;  therefore,  to  minimize  this  effect,  t  must  be  made  as 
small  as  possible.  Also,  the  output  will  respond  to  within  5  percent 
error  in  a  delay  time,  Tj),  of  approximately  r. 


SECTION  V 
RAMP  INPUT 


In  order  to  theoretically  analyze  the  error  of  the  V/F  converter 
data  acquisition  system  in  the  data  recovery  of  a  ramp  input,  time, 
t  =  0,  is  set  at  the  start  of  the  summing  interval  in  which  the  ramp 
function  begins. 

The  ramp  function  is  defined  as : 

l<“  V 

fi  (t  —  a)  =  <  (21) 

lt-a  l>  a 

For  a  normalized  ramp  function  which  has  a  rise  time  of  /3  summing 
intervals,  the  input  function  is  expressed  as: 


Vi(t)  -  J-  [/I,  (t  -  Kr)  -  /i2  (t  -  (K  P)  r)] 
Pr 


(22) 


The  value  of  VTq  is  dependent  upon  the  value  of  K  +  /3.  For 
(K  +  /3)  <  1,  the  combination  of  Eqs.  (3)  and  (22)  becomes: 

Vro  =  —1—  r  [ft,  (t  -  Kr)  -  pi,  (t  -  (K  +  p)  r)]  dt 


i  r{K^)T  i  rr 

=  I  (t  -  Kr)  dt  +  -  I  dt 

^  -Kr  *|C  K  +  P^t 

=  '  ['At2  -  Kn](kT+>-[t]r 

Pt1  Kr  :T  ( K  +  /S)r 


=  _L  ['A  (K  py  7*  -  'AK1  T2  -  K  (K  +  p)  7*  +  K*  r2] 

P*1 

+  1  [r  -  (K  +  j8)  r] 

=  L  +  1  -  (K  -  p) 

-  1  -  K  -  A  <23) 

for  0  <  K  +  0  <  1. 
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For  (K  +  0)  >  1,  VTq  is  defined  by  the  integral: 

v,‘  '  wS'il ' Hd* 

=  _2_[ 'At2  -  Krt]'  =  -i-  [H  -  KK1  -  K.  +  K2] 
j Sr1  /Sr2 

(1  -  K)2  (24) 

2/S 

for  (K  +  /3)  >  1. 


For  all  summing  intervals  up  to,  but  not  including,  the  one  in  which 
the  ramp  breaks  to  a  zero  slope,  the  value  of  the  output  is  equal  to  the 
value  of  the  input  at  the  midpoint  of  the  summing  interval.  That  is: 


for  0  <  n  <  Int  (K  +  0). 


w„  .  v,  (*±J.  ,) 


For  the  interval  n  =  Int  (K  +  |3),  the  output  value  is: 

J(K  +j8)r  (n-l)r 

(t-Kr)dt  +  lf  dt 
nr  +j8)r 


jsV  J  (t  ~  Kf)  dt  +  r  f 

nr  j(K+8)t 


]  r„  (K+/8)r  ,  (n  +  l)r 

~-VAl2  -  Krt]  +  I  [t] 


<K+jS)r 


=  !S(K  +  8)2  -  ~  +  Jjj-  +(n  -  §  -  (K  +  /3) 

=  (n  +.  ])  -  K  -  8/2  - 

V  2/S 

Figure  10  depicts  the  input-output  relationships  for  a  ramp 
function. 


-i  OK  r£  Jr  5i-  t 

2  2  2  2 
Fig.  10  Input-Output  Waveforms  for  Ramp  Input 


(25) 
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5.1  LINEAR  INTERPOLATION  ERROR  ANALYSIS  OF  RAMP  INPUT 

As  shown  in  Fig.  10,  the  emax  occurs  at  either  t  =  Kr  or  t  =  (K  +  |3)t; 
therefore,  emax  is  equal  to  the  maximum  of  these  two  values.  An  expres¬ 
sion  for  the  emax  is  so  cumbersome  to  deal  with  manually  that  a  digital 
computer  program,  depicted  in  Table  I-IV  (Appendix),  was  written  to 
calculate  the  many  solutions  required  to  determine  the  minimum  j3  for 
5  percent  or  less  error.  The  program  calculates  the  0min  for 
emax  -  ^  percent;  then  j3  is  varied  from  0. 1  to  j3min  in  increments  of 
0.  1  for  K  =  0.  0  and  0.  50  and  the  value  of  emax  determined  for  each  set 
of  parameters.  The  family  of  curves  that  was  calculated  is  depicted  in 
Fig.  11.  These  two  values  of  K  are  the  extreme  cases;  all  other  values 
are  in  between.  The  plots  in  Fig.  11  indicate  that  a  value  of  j3  >  5.  0 
assures  a  response  within  5  percent. 


Fig.  1 1  Plot  of  emax  versus  j3  for  Ramp  Input  with  K  Constant 


5.2  DATA  POINT  ERROR  ANALYSIS  OF  RAMP  INPUT 


The  value  of  Emax,  as  shown  in  Fig.  10,  occurs  at  either  ti  =  l/2r 
or  t2  =  [Int  (K  +  /3)  +  1/2]  t.  The  expression  for  Emax  is  the  maximum  of: 


E*  =  100  | Vj  (Hr)  -  VrJ 
or 

E2  =  100  l\'i(^±-!r)  -  Vr„' 

where  n  =  Int  (K  +  /3). 
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Fig.  12  Plot  of  Em<JX  versos  jSfor  Ramp  Input  with  K  Constant 

The  digital  computer  program  that  calculates  emax  was  written  for 
the  Emax  option  also.  The  value  of  0min  for  Emax  <  5  percent  is  first 
calculated;  then  is  varied  from  0.  1  to  $niin  in  increments  of  0.  1  for 
K  =  0.  0  and  0.  5,  since  these  errors  are  the  extreme  cases,  and  the 
value  of  Emax  determined  for  each  set  of  parameters.  Figure  12  depicts 
the  family  of  curves  that  was  calculated.  These  plots  indicate  that,  for 
a  value  of  £3  >  2.  5,  the  data  points  are  in  error  by  5  percent  or  less.  To 
illustrate,  for  a  summing  interval  of  2  msec,  any  ramp  with  a  rise  time 
of  5  msec  or  greater  will  be  reconstructed  with  all  data  points  within 
5  percent  error. 


SECTION  VI 

SUMMARY  AND  CONCLUSIONS 


6.1  SUMMARY 

A  voltage -to -frequency  data  acquisition  system  was  theoretically 
analyzed  for  errors  that  occur  when  it  is  used  as  a  sampling  system. 
Since  the  standard  sampling  theorem  is  not  directly  applicable,  the 
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system  was  simulated  on  a  digital  computer  to  aid  in  the  determination 
of  parameter  criteria  for  expected  levels  of  accuracy  for  sinusoidal, 
step,  and  ramp  inputs. 

Parameter  limits  were  determined  for  such  parameters  as  m, 
summing  intervals  per  period  of  a  sinusoid,  K,  delay  of  step  and  ramp 
inputs  with  respect  to  start  of  summing  interval,  and  0,  ratio  of  ramp 
rise  time  to  summing  interval  length. 

An  investigation  of  the  effect  of  signal  synchronization  was  accom¬ 
plished  for  the  sinusoid,  and  an  experiment  was  performed,  using  a 
sinusoid,  to  determine  the  validity  of  the  theoretical  solutions. 


6.2  CONCLUSIONS 

The  voltage -to-frequency  conversion  for  data  acquisition  does  not 
offer  accurate  recovery  (5  percent  error)  of  sinusoidal  data  until  the 
ratio  of  summing  interval  to  period  is  about  11.4,  i.  e. ,  11.4  summing 
intervals  per  period.  This  means  that  for  a  2 -msec  summing  interval, 
any  frequency  component  above  43.  9  Hz  has  a  maximum  error  greater 
than  5  percent  and,  from  Fig.  3,  any  component  of  151.  5  Hz  or  greater 
is  in  error  at  its  maximum  by  more  than  50  percent. 

The  voltage-to-frequency  conversion  system  has  an  inherent  error 
because  all  data  points  do  not  coincide  with  the  input  waveform.  For 
the  sinusoid,  the  ratio  of  summing  intervals  to  period  must  be  12.  8  for 
the  data  point  to  be  within  1  percent  of  the  input  waveform.  At  the  value 
of  m  =  11.4,  emax  was  4.  95  percent  and  Emax  is  1.  26  percent;  there¬ 
fore,  a  sizeable  portion  of  the  error  is  the  result  of  the  data  point  error. 

The  step  input  error  is  independent  of  the  summing  interval  and  is 
symmetrical  with  respect  to  K  =0.5  for  both  emax  and  Emax.  However, 
delay  time,  time  required  to  bring  the  output  to  be  within  5  percent  of  . 
the  input,  is  a  function  of  K  and  r,  and  as  K  varies,  the  time  delay 
varies  in  the  range  from  0.  45t  to  1.  05t. 

Ramp  inputs  can  be  followed  within  5  percent  if  the  ratio  of  ramp 
rise  time  to  summing  interval  is  larger  than  five.  This  is  to  say  that 
a  2 -msec  summing  interval  can  accurately  follow  a  ramp  of  rise  time 
greater  than  10  msec. 

The  ramp  function  has  an  unusual  Emax  versus  0  for  constant  K 
relationship.  At  low  and  high  values  of  K,  the  Emax  oscillates  as  a 
function  of  |3.  At  midvalues  of  K,  the  relationship  is  nearly  exponential. 
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If  j3  >  12.  5,  all  data  points  are  within  1  percent  of  the  input  waveform. 
The  reverse  conditions  exist  in  the  relationship  of  emax  versus  $  for 
constant  K. 

An  additional  error  that  may  be  significant  at  small  summing 
interval  lengths  is  resolution  error  from  the  low  pulse  counts  that 
occur.  However,  errors  from  resolution  were  not  included  in  this 
analysis. 
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APPENDIX 

COMPUTER  PROGRAMS 
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TABLE  l-l 

COMPUTER  PROGRAM  TO  FIND  Emox  AND  emox  FOR  A  SINUSOIDAL  INPUT 
AS  m  VARIES  FROM  1.0  TO  25.0  IN  INCREMENTS  OF  0.1 


L. 00 01  /JOB  GO, T IME=10  ”7‘-  ■  ~  '  - 

L. 0002  DIMENSION  XI-1A<  241),DI  FMAC  241 )  -  •-  --  .  ,■ 

L. 0003  C  I ND=0  (E(MAX)  FOR  DATA  POINT)  .  .  '  > 

L. 0004  C  I ND=1  (E(MAX)  FOR  LINEAR  INTERPOLATION  AT  ANY  POINT) 

L. 0005  READC  5,1)  1  NO 

L.  0006  1  FORIIAT(IIO) 

L. 0007  PSI=0.  -  -  v  . 

L . 0008  DO  103  IM*10,  250  ,  V  '  ^ ',r  >• 

L. 0009  I  =111-9  -  .  •  ,  1 

L.  0010  XM=FLOAT(  IID/10. 

L.0011  XIIA(I)-XI1 

L.0012  DIFMACI  )=VEI1AX(  XII,  PS  I  ,  IND) 

L.0013  103  CONTINUE.  - 

L.  0014  WR I  TE(  6, 2 )  '  .  •  ,  • 

L.  0015  2  FORMATCl'-)  .  '  ,  '  .  *  :  -V 

L. 0016  I F  C I ND)100, 100,101 

L. 0017  100  WRITE(6,3) 

L.0018  3  F0UI1ATC lOX'DATA  POINT  ERROR  ANALYSIS'/) 

L. 0019  101  WRITE(G,4)  { XMAC  J  ),D I  FMAC  J ), J =1,241) 

L.0020  4  FORMAT (4{5«,M,4XI  ECMAX),1  4X)//4(2F8.2,4X))‘  ' 

L.0021  STOP,  -  .  1  *  .  i 

L. 0022  END 
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TABLE  l-l  (Concluded) 


L.  0001 

FUNCTION  VEHAX(XH,PSI,IND) 

L. 0002 

DIMENSION  FAC  252) ,TA C  25  2 ) 

L.  0003 

PI =  3 . 1415926 

L. 0004 

P5I =PSI *PI/180. 

L. 0005 

IFCXM-25. )99,99,93 

L.  0006 

98 

STOP  2525 

L.  0007 

99 

DO  110  1=10,100,10 

L. 0008 

XP=0.QQ1 

L. 0009 

DO  109  1=1,2  . 

L. 0010 

1 XC=XM*F  LOAT ( L)+XP 

L.  0011 

1 F  C IXC-C IXC/ 10) *10) 103, 103, 110 

L . 0012 

103 

XP=XP+0 . 9 

L. 0013 

109 

CONTINUE 

L. 0014 

GO  TO  111  -  . 

L.0015 

110 

CONTINUE 

L. 0016 

111 

NS=( 1 XC+25)/10 

L.  0017 

ARG1=PI/XM 

L-.0018 

EMAX=0. 

L.  0019 

DO  100  N=l, NS 

L. 0020 

XN-N-2 

L. 0021 

ALPIIA=2.*XN+1.  :  - 

L.0022 

TAC  N  )  =ALPI  IA/XI1/2  . 

L.0023 

ARG2=ALPIIA*ARG1+PSI 

L .  0024 

FAC  N  )  =S 1 NC ARG2 ) *S 1 NC ARG1 )/ARGl 

L.  0025 

ERR  =A8S { S 1 N  C  ARG2 ) -FA  C  M ) ) 

L.0026 

EMAX=AMAX1( ERR,EMAX) 

L.0027 

100 

CONTINUE 

L.0028 

IFC 1 ND)103,1Q3,101 

L.  0029 

101 

EIIAX=0. 

L.0030 

IXS=TACNS)*100.+1. 

L  -  0031 

DO  102  1 X-l, 1 XS 

L.  0032 

X=F  LOAT (IX)/ 100 . 

L.0033 

TEST=TAB2CFA,TA,X,NS) 

L. 0034 

DEV=A0SC  SI HC  2. *PI *X+P5 1 )-TE5T) 

L. 0035 

OIAX=AMAXl(EMAX,DEV) 

L. 0036 

102 

CONTINUE 

L.  0037 

103 

VEMAX  =  EMAX*100 . 

L.  0038 

RETURN  •  1  *  ' 

L*  0039 

END. 
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TABLE  Ml 

COMPUTER  PROGRAM  TO  FIND  THE  LARGEST,  SMALLEST,  AND  DELTA 
Emax  OR  emax  as  ij,  VARIES  FROM  0  TO  180  deg  IN  INCREMENTS  OF 
1  deg  AND  m  VARIES  FROM  1  TO  25  IN  INCREMENTS  OF  0.5 


0001 

_  0.002 _ 

0003 

JQQ-0A. _ 

0005 

.0006 _ 

0007 

0008 

0009 

J3Q.1.Q. _ 

0011 

0012 

0013 _ 

00  14 

.00.15 _ 

0016 

0017.  . 

0018 

0019  _  _ 

0020 

.0021 _ 

0022 
0023 
CO  24 
0025 
0026 

.0027.  _ 

0028 
0029 
00  30 
00  31  . 

0032 

C033  _ 

0034 

00  35.  _ 

0036 

0037  .  . 

00  38 

0039  . 

0040 

0041 _ 

0042 

0043 

0044 

0045_ 

0046 

.  .Q.Q&L. .  . 
0048 
0049 

0050 


DIMENSION  TA1252),FA(252) 

_ P.1 53.  ..1.4 1.59 2  6 _ _ _ 

DO  906  ND=1 • 2 

_ JW  0=N.QrJ _ 

WRITE (6,1) 

_  L_  FORMAT  ( ♦  1 »  ) _ _ 

IF! INDJ902, 902,903 

902_WRITE(6,_2) _ _ _ _ _ _ ... _ 

2  FORMAT {//lSX'OATA  POINT  ERROR  ANALYSIS*) 

903_WR_IiE.L6x.3J _ 

3  FORMAT I/23X*LARGEST* 13X* SMALLEST  DELTA* /7X 

_  1  *  M ♦ 2 X2J  6X*PSI*5X»E(MAX)«) »4X*  E  t  MAX) ♦) _ 

DO  906  IM=10,250,5 

_ XM=FLJ3  Aid M  )  / 1 0 . _ 

DO  99  1  =  1,252 

99  TAtH=l2.*FL0AT(  I-2)  +  l.)/XH/2. 

EMAOUT=0. 

_  EMIOUT=L50_._ _ _ 

00  102  L=l6,100,10 

_  XP=0_.001 _ . _ _ 

DO  i'dl  1  =  1,2 

i XC=XM*FLQATIL)+XP _ 

IF  Cl XC-t IXC/10)*10) 100, 100, 102 

100  XP=XP+G.9 _ . _ 

101  CONTINUE 

GO  TO  103  _ 

102  CONTINUE 

JO 3  NS^CIXC  +  251/iQ _ 

ARG 1 =P I / XM 

DO  900  1  =  1,181 _ _ _ 

CSl=I-i 

_PSI=CSI*PJ/J80. _ _ 

EMAX=0 • 

_ DO  104  _N=J  iNS _ 

ARG2=l 2 . * FLOAT ( N-2 ) + 1 . ) *ARG 1*PS I 

_  FA(N)  =  S1N( ARG2)*SIN(ARG1)/ARG1 _ 

ERR=ABSi  SIN  I  ARG2I-FAIN)  )~*100. 
EMAX=AMAX1(ERR,EMAX) _ 

104  CONTINUE 

IF(  IND)  107,  107,105 _ 

105  emax=o. 

IXS=TAINS_)*100_. _ 

K=2 

_  00  106  IX=_1,_IXS _ 

X=FLUAT ( IX-i)/100. 

70J  PI F=X-TA I K I _ 

IFIDIF) 703, 703,702 

702  K=K+1 _ 

GO  TO  701 

703  DEV=ABSISIN(2.*PI*X+PSI )-DIF*(FAlK)-FAIK-l) ) 
1/(TA(K)-TA{K-1 ) j-FAIK) 1*100. 

EM AX=AM AX1 1  DEV , EMAX ) 
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TABLE  l-ll  (Concluded) 


CO  51 

106 

CONTINUE 

0052 

107 

IF(EMACUT-EMAXI8CO,0O1,3O1 

OC  53 

000 

EMA.1UT  =  EMAX 

CC  54 

PS IL=CS I 

OC  55 

001 

IF(CMAX-EMIUUT)802* 9G0 , 90C 

00  56 

00  2 

EMinUT=EHAX  .  . 

00  57 

PSIS=CSI 

0050 

900 

CONTINUE 

0050 

DIFM4=EMA0UT-EMI0UT 

CC60 

WRITE!  6, 4)  XM,PSIt-»EMAOUT»PSlS»EMIOUT,DlFMA 

0061 

4 

FORMAT (6F10.  ?) 

0062 

906 

CONTINUE  .  .  _ 

OC  6  3 

STOP 

0064 

END 
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TABLE  Mil 

COMPUTER  PROGRAM  TO. COMPARE  THEORETICAL  AND  ACTUAL 
•mox  or  Emax  FOR  SINUSOIDAL  INPUT 


L.0001 

/JOfi 

fil),TlliE=10 

L. 0002 

1)1  MEL'S  1011  UATAC292),FAC292),TAC292),FCC140),T8C140), 

L.000 3 

1XMC  12)  ,  El 1AX (  12 >  #  ET.-IAX ( 12  ) 

L. 0004 

11 EAD (5,1)  1  IK),  IZC,  IUC,VU,-)A7A 

L.0005 

1 

rOUMAT{3U0,riO.3/ClGr5.1>) 

L. 0000 

UUITE(G,2) 

L.0007 

2 

FORIIATC 1 1  * ) 

L.000S 

TAU =0.002 

L.OOO'J 

VKT=TAU*F  LOAT(  1 UC-  1  ZO/VU 

L. 0010 

ZF=FL0AT( 1 ZC)*TAU 

L.0011 

vma::=o. 

L .0012 

1)0  100  11=1,292 

L. 0013 

TA  ( 11 )  =F  LOATC  2*11-  2  7 )  *TAU/  2 . 

L.0014 

facii)=cuatach)-zf)/vi;t 

L. 0015 

10U 

viia:c-amaxi(v;iax,ass(fa(:o)> 

L. U01G 

1)0  104  0  =  2,13 

L. 0017 

XJ=15-J 

L. 0012 

TAUJ  =TAU*X J 

L.OOl'J 

VICTJ=VICT*XJ 

L. 0020 

TZCJ=ZF*XJ 

L. 0021 

X!  1  ( J  - 1 )  =0 . 0  5  /  TAU  J 

L . 0022 

et:iax(j-i)=vemax(xiko-i),o.,  i ul> > 

L. 0023 

EUR  =0 . 

L.O024 

11=279/ ( 15-J )+l 

L.0025 

DO  102  .L-1,K 

L. 0020 

UM=0. 

L. 0027 

LL=l4+(15-0)*(L-2) 

L. 0028 

LU=LL+14-d 

L. 0029 

1)0  101  1  =LL,  LU 

L.U030 

101 

Oll=UI!+l)ATA(  1  ) 

L. 0031 

fb(l)=(oh-tzco)/vi;tj 

L.0032 

TBCL)=FLOATC2*L-3)*TAUJ/2. 

L.0033 

VV1»7AB2( FA,  TA,  TR(  L) , 292 ) 

L. 0054 

r.  11 11  =AI  1AX 1  ( AI5S  ( V  V 1-  F!»  (•!_)),  ERR ) 

L.  01)35 

102 

COI1T1  HUE 

L.003G 

1  F  C  1110)105,  105,  10G 

L .  0037 

10G 

ERR=0. 

L.0038 

DO  103  1=1,1001 

L . 0030 

TIME=F LOATC  l-l)*i). 0005 

L. 0040 

VV1=7A32C  FA, TA, TIME,  292 ) 

L.0041 

VV2-TAB2(FC,TB,TIME,IO 

L. 0042 

103 

EUR=AMAX1CAB5( VV1-VV2) ,  ERR) 

L .0043 

105 

EiMAX  ( J  - 1  >  =  EUR/  VMAX  *  1 0  0 . 

L . 0044 

104 

COI1TIIJUE 

L, 0045 

IFC  1110)107, 107, 100 

L.004G 

107 

Will  TEC  G,  3) 

L.  0047 

3 

FORMAT Cl OX1 DATA  POIUT  ERROR  ANALYSIS’) 

L.0048 

10Q 

Will  TEC  0,4)  C  XM(  J ),  EMAXC  J ),  ETI1AXC  J  ),  J=l,  12) 

L. 0049 

4 

FORMATC 2 C20X’ ACTUAL' 2X ' THEORY ' 4X )/ 2 C 15X 'M ' 4X ' ECMAX) ' 2X 

L, 0050 

1 ' EC MAX) ,4X)//2C10X3F3.2,4X)) 

L.0051 

STOP 

L. 0052 

EIID 
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TABLE  MV 

COMPUTER  PROGRAM  TO,  FIND  ema^OR  Emox  FOR  A  CONSTANT  K  AS  0 
VARIES  FROM  0.1  TO  12.0  IN  INCREMENTS  OF  0.1 


L.0001 

/JOB 

no 

L . 0002 

n II1RNS 1  OH  F.RR(  150),.B.(15P) . 

. 

L.0G03 

C 

1  00=0  (  E(MAX)  FOR  DATA  PO 1  NT) 

L. 0004 

C 

IND  =  1  ( E(MAX)  FOR  LI  HEAR  INTERPOLATION  AT 

ANY 

POINT) 

L. 00  05 

DC)  105  0  =  1.2 

L.  00 06 

1 ND=2- J 

L.  0007 

Will  TE£  6, 2) 

L. 0008 

2 

FORMAT ( '  1' ) . . . .  _  ... 

.... 

_ _  .. 

L. 0009 

IF(  1110)97,1)7,38 

L.0010 

97 

1  P=21013fi2752 

L.  0011 

00  TO  99 

L. 0012 

98 

1  P=107795257G  .  . 

L. 0013 

99 

NB=120 

L.0014 

nn  106  1=1^4 

L.  0015 

XK-FI.OATC  l-D/4. 

L.  0016 

DO  104  10=1.1(0  .... 

_ _ _ _ _ 

L.0017 

B(  1 B) =F  LOAT ( IB)*0. 1 

L.  00 18 

ERR(  IB)=EI1AXR(B(  IBJ.XK.  IND) 

. 

_  .  . 

L.  00 19 

104 

CONTI  IHJE 

L.  002  0 

UR  1  TE(  6.3)  1  P,  XK.  ( B (.1  fllv-ER3t.l-B.)..-l.B=-l,-NB.) .... 

_ _  _ 

..  .  _. 

L.  0021 

3 

FORMAT ( // '  ERROR  ANALYSIS  FOR  BETA*  Al*  KITH  K  = 

>  'F  5. 2 

L. 0022 

l//6( 4X' BETA  E ( MAX ) 1 4X) // 6 ( 2F  8 . 2. 4X ) ) - ‘ 

L.  0023 

10  5 

CONTINUE 

L. 0024 

STOP 

L.  002  5 

END 

L.  0026 

FUNCTION  EMAXR ( B ET A..XK.  I  NR )..  .  . 

■  ...  _ _ 

_ _ _ _ 

L.0027 

DIMENSION  FIK  23)  ,TfM  23) 

L. 0028 

C 

1 ND=0  ( E( MAX)  FOR  DATA  POINT) 

L .  0029 

C 

1  NO  =  1  (F.(MAX)  FOR  LINEAR  INTERPOLATION  AT 

AMY 

POIMT) 

L. 0030 

RAMPF  ( T)  =AMI  f.‘l(  1. .  AflAXl(T-XK,  0, )  /BETA) 

L . 0031 

FN( 1) =0. 0 

L. 0032 

TNC  1  )=-!).  5 

L.0033 

DO  100  J  =  2 . 2 3 

L . 0034 

TIM  J)=TN(  J-D  +  15  .... 

... 

_ 

L.0035 

100 

FN( J) =1. 0 

L. 0036 

TS1=XK  +  BET.A  ..  . 

_ 

. 

L. 0037 

IS1=TS1 

L.0033 

TS  2  =  1 S 1 

L. 0039 

IS2-IS1+1 

L.0040 

TS3  =TS2-XK 

L. 0041 

FN(2)=1.-XK- BETA/2. 

- 

L. 0042 

IF( 131-1)105.101,101 

L. 0043 

101 

FN( 2) =( l.-XK)**2/BETA/2. 

L. 0044 

IFf IS1-711 04 ,102,102 

L. 0045 

102 

DO  103  J=3, 1 S2 

L. 0046 

103 

FIM  J)  =RAHPF  ( TN(  >J) ) 

.  . 

L.0047 

104 

FN(  IS2+1)  "TS3+1.  -(BF.TA+TSj**2/RETA)/2. 

L. 0048 

105 

IFC IND)106, 106, 107 

L. 0049 

106 

EMl=ABS(RAMPF(0.5)-r;.(2)) 

L . 0050 

FM2  =ABS  C  RAMPF  CTS2+0. 5) -FIM  IS2  +  1)) 

L. 0051 

00  TO  108 

L.0052 

10  7 

EI11=ABS  (  TAB2  ( FN,  Til,  XK..13J  )  .  .  ..  .  .. 

L. 0053 

EII2=ABS ( TAR 2 ( FN.TN,  TS1,  13)  - 1 . ) 

L. 0054 

108 

EMAXR  =AI  SAX  1  (  EMI,  EM2  )  *10  0 . 

L. 0055 

RETURN 

1.0056 

END 
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